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Materials adsorbed to the surface of a fluid - for instance, crude oil, biogenic slicks, or 
industrial/medical surfactants - will move in response to surface waves. Due to the diffi¬ 
culty of non-invasive measurement of the spatial distribution of a molecular monolayer, 
little is known about the dynamics that couple the surface waves and the evolving den¬ 
sity field. Here, we report measurements of the spatiotemporal dynamics of the density 
held of an insoluble surfactant driven by gravity-capillary waves in a shallow cylindri¬ 
cal container. Standing Faraday waves and traveling waves generated by the meniscus 
are superimposed to create a non-trivial surfactant density field. We measure both the 
height held of the surface using moire-imaging, and the density held of the surfactant via 
the huorescence of NBD-tagged phosphatidylcholine, a lipid. Through phase-averaging 
stroboscopically-acquired images of the density held, we determine that the surfactant 
accumulates on the leading edge of the traveling meniscus waves and in the troughs of 
the standing Faraday waves. We ht the spatiotemporal variations in the two helds using 
an ansatz consisting of a superposition of Bessel functions, and report measurements of 
the wavenumbers and energy damping factors associated with the meniscus and Fara¬ 
day waves, as well as the spatial and temporal phase shifts between them. While these 
measurements are largely consistent for both types of waves and both helds, it is notable 
that the damping factors for height and surfactant in the meniscus waves do not agree, 
suggesting the presence of longitudinal waves. 


1. Introduction 

The calming effect of surface oil on oceanic waves was recognised by sponge and pearl 


divers as early as the time of Pliny the Elder. Franklin et al. (17741, citing Pliny’s ac¬ 
count, also reported on sailors calming ocean waves with oil. Franklin concluded from 
experiments that oil significantly damps short-wavelength, small-amplitude waves, and 
conjectured that it would also induce damping of larger long wavelength waves. Lord 
Kelvin (Thomson 1871) hrst identihed short-wavelength capillary-driven waves and long- 
wavelength gravity-driven waves as two regimes within the same mathematical formu¬ 
lation. Kelvin’s famous dispersion relation describing the wavelength-dependence of the 
speed of gravity-capillary waves shows a crossover from capillary to gravity waves at the 
minimum wave speed. 

Much effort has since been devoted to studying the damping effect of surfactants on 
gravity-capillary waves (|Reynolds 1880 Levich 1941 Dorrestein 1951 Goodrich 1961 


Lucassen-Reynders & Lucassen 1970|). With or without a surfactant, energy is lost due 


to the viscosity of the fluid and the vorticity of the flow. By adding a surfactant, the 
surface can experience additional tangential stresses which, balanced by viscous stress in 
the bulk, increase the vorticity of the bulk flow. These tangential surface stresses are a 
consequence of the complex relation between surface stress and strain rate, which depends 
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upon the bulk flow, surface velocity, and surfactant density. The primary contribution 
to the tangential surface stress is the Marangoni effect (Behroozi et al. 2007), in which 
the advection of the surfactant monolayer by the passing wave results in gradients in the 
surfactant density field, inducing surface tension gradients and vorticity in the bulk fluid 
flow. Thus, beyond quantifying the damping effect of surfactants on gravity-capillary 
waves, it is important to also understand the distribution of surfactant on the surface 
of the fluid. In this paper, we combine two experimental techniques to measure both 
the surface height and surfactant density fields. Using a physically-motivated ansatz, we 
quantify both fields and measure parameters such as the complex wavenumber and phase. 
This provides a full spatiotemporal description, including the locations of surfactant 
accumulation within the wave pattern. 

Surface waves can be generated and sustained in a variety of ways, depending on 
how energy is injected into the fluid system. A vertically-vibrated fluid will generate 
both traveling meniscus waves, excited by the contact line at the container wall, and 
standing Faraday waves, first observed by Faraday (1831). When the driving acceleration 
a(t) = ao sin (cot) is weak, only meniscus waves perturb the fluid surface. However, when 
do is increased beyond a critical acceleration amplitude a c , so that the injection of energy 
exceeds the dissipation due to the bulk viscosity, the meniscus wave becomes unstable to 


Faraday waves which grow to a finite amplitude (Benjamin & Ursell 1954 Chen & Vinals 


1999). The parameter a c depends on the driving frequency w, the container geometry, 
and fluid properties (Douady|1990 Edwards fe Fauve|l994 Bechhoefer et al. 1995| Chen 


& Vinals 1999). The emergent Faraday waves can be either harmonic (frequency w) or 


subharmonic (frequency to/ 2), while the meniscus waves are always harmonic. 

Recently, there have been advances in our understanding of the relationship be¬ 
tween the surface height field h( r, t) of meniscus waves and its interaction with sur¬ 
factant molecules adsorbed to the surface, described by the density field T(r,t). The 
theory of surfactant-laden meniscus waves in a cylindrical geometry, developed by |Bock 


(1991), Saylor et al. (20001, and Picard & Davoust (2006), proceeds from the treatment 


of 2-dimensional traveling gravity-capillary waves in a Cartesian geometry (Lucassen- 
|Reynders fe Lucassen|[l970 1. In this theory, the fluid motion in the incompressible bulk 
is modeled with the linearized Navier-Stokes equations, and the vertical and horizontal 
displacements of the surface satisfy the surface stress and kinematic boundary condi¬ 
tions. Bock (1991) showed that for a fluid of equilibrium height ho, the deviation in the 
surface height field Ah(r, t ) = h( r, t ) — h 0 for inward and outward traveling cylindrical 
gravity-capillary waves follow Hankel modes of the first and second kinds. Recognizing 


that meniscus waves are a superposition of inward and outward traveling waves, Saylor 


et al. (2000) found that meniscus waves follow a Jq Bessel mode. Subsequently, Picard 


& Davoust (2006) derived an expression for the deviation of the surfactant density field 


from an equilibrium density r 0 . The quantity Ar(r,f) = T(r, t) — r 0 also follows a J 0 
Bessel mode, with the same wavenumber and damping factor as the surface height field. 
Because the motion of a fluid surface element depends upon the surface compression 
modulus e = Toyf, the magnitude of e controls the phase shift between Ah and AT 
fields (Lucassen-Reynders & Lucassen 1970). 

For Faraday waves, Kumar & Matar (J2002, 2004) used a linear stability analysis to 
determine the wave’s critical acceleration and wavelength and predicted that the Faraday 
wave emerges with a spatial displacement (phase shift) between the fields Ah and AT. 
Inspired by reports (Douady et al. 1989|) of a rotating Faraday wave, Martin & Vega 


(2006) conjectured that the rotation was induced by the presence of a contaminant and 
showed that the drift can be generated when the reflection symmetry of the streaming 
flow breaks. Ubal et al. (2005) executed numerical simulations of the gravity-modulated 
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Figure 1. (colour online) Experimental apparatus. The container (black) is connected by 
low-friction air bushings to an electromagnetic shaker, and filled with a thin layer of water 
onto which a monolayer of NBD-PC surfactant is deposited. The shaker vibrates the container 
vertically with acceleration a(t) = ero sintaf, generating meniscus and Faraday waves. The flu¬ 
orescence imaging system involves a ring of blue LEDs around the outer edge that excites the 
fluorophore on the lipids; digital images of the fluorescence emission are collected through an 
optical notch filter. A silicon wafer sits flush with the bottom surface of the container, and re¬ 
flects the excess light from the LEDs away from the camera. The moire imaging system consists 
of a patterned light source shining downward on the fluid surface, and the camera which records 
the distorted reflection of the pattern. 


Navier-Stokes equations and predicted that neither Ah nor AT evolve sinusoidally in 
time. They found that Ah lags behind AT, characterised as a temporal phase shift 
between the two fields. 

Many experiments with surfactant-laden gravity-capillary waves have focused on quan¬ 
tifying energy dissipation through either direct or indirect measurements of the surface 
height field alone (Case fe Parkinson||1957 Davies & Vosc 1965; Jiang et al. 1993 Hen¬ 


derson & Miles 1994 Henderson 1998 Saylor et al. 2000; Behroozi et al. 2007), and 
these works have demonstrated that the damping rate does not depend on whether the 
gravity-capillary waves are either standing or traveling. Despite this progress, a quanti¬ 
tative understanding of the spatiotemporal dynamics of T(r, t) remains elusive. Previous 
measurements of T have utilized surface potential measurements in a small region near 
an electric probe (Lange & Hiihnerfuss 1984 Hiihnerfuss et a7|l9 85; Lange & Hiihncrfuss 
1986). These studies report a temporal phase shift between AT and Ah, with the caveat 


that it is necessary to account for both the response time of the probes and the finite 
time for the wave to pass from the surface potential probe to the wave height probe 


Lucassen-Reynders (1987). In this paper, we present novel methods for measuring the 


spatial and temporal behaviors of AT and Ah, using stroboscopic imaging to avoid this 
complication. These methods allow us to quantify the accumulation of surfactant rela¬ 
tive to the crests/trough of the surface waves. In addition, we test the extent to which 
the spatiotemporal dynamics can be described by a superposition of solutions for the 
meniscus and Faraday waves. 

Our experimental apparatus, shown schematically in figure |T] and described in detail 
in © consists of a shallow cylindrical container holding a thin layer of water covered 
with a monolayer of fluorescently-tagged lipids. We drive the system with vertical oscil¬ 
lations just above onset for the Faraday instability and make stroboscopic measurements 
of both the surface height field h(r,t) and the surfactant density field r(r,f). To obtain 
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h , we illuminate the fluid with a target pattern, and use a combination of ray-tracing 
and nonlinear fitting to invert the resulting moire image. For T, we phase-average the 
fluorescence intensity from the tagged molecules using images acquired stroboscopically 
over many cycles of the driving oscillation, and convert these to quantitative measure¬ 
ments of the surfactant density field. A key advantage of these techniques is that they are 
non-invasive. Additionally, this and similar optical techniques ( Vogel et ai. ||2001 Fallest 
et aL|20i0 Strickland et a/.|2014 Swanson et al.\‘2(y\A ) for measuring T depend upon the 
mean distance between the fluorophores and therefore are not affected by the dynamics 
of either molecular rearrangement or domain formation/relaxation. 

In © we show that our data are well-approximated by a linear superposition of a 
standing Faraday wave mode and a traveling meniscus wave mode. We decompose the 
data into these two separate components in order to examine the spatiotemporal dy¬ 
namics of each. We determine that the surfactant accumulates on the leading edge of 
the meniscus waves and in the troughs of the Faraday waves. The meniscus waves are 
represented by Jq Bessel functions. The fields Ah and AT have the same wavenumber 
within experimental error, but the damping factors for the two fields do not agree. Since 
the temporal phase shift between the two fields is measured to be as large as 2 radians, 
in excess of the theoretical maximum of n/2 dLucassen-Reynders & Lucassenjl970j), we 
conjecture that longitudinal waves (Lucassen 19686) may also be present. The Faraday 
waves are represented by J„ Bessel functions. The fields Ah and AT have the same sym¬ 
metry number n and the same wavenumber within experimental error. We observe that 
the whole Faraday wave pattern rotates around the center of the pattern in the same 
direction as an observed spatial phase shift between Ah and AT fields. Both fields evolve 
sinusoidally and are temporally phase shifted by roughly 2.4 radians. 

In (|4j we discuss our observations about surfactant-covered meniscus and Faraday 
waves in context of their respective theoretical frameworks. We also consider the possible 
presence of resonantly-excited longitudinal waves in our system, which could account for 
the anomalously large temporal phase shifts and the disagreement in the damping factors. 
We highlight the suitability of using measurements of AT and the temporal phase shift 
between Ah and AT as a way to probe the interfacial rheology. Finally, we contrast the 
dynamics of these molecular monolayers with the dynamics of a monolayer of millimetric 
sized particles (Sanl~eFaT||2014). 


2. Experiment 

We excite traveling meniscus waves and standing Faraday waves on a surfactant- 
covered fluid layer by subjecting the system to a vertical sinusoidal oscillation, shown 
schematically in figure [l] A mechanical driving system provides the vertical acceleration, 
and two quasi-independent imaging systems measure the response of both the surfactant 
and the fluid surface. The fluorescence imaging (FI) system measures the surfactant den¬ 
sity field r(r,f) (details provided in { 2.21, and the moire imaging (MI) system measures 
the surface height field (j |2.3[ ). Finally, we use a physically-motivated ansatz ((2.4) to 
support the ray-tracing and image analysis in order to convert the MI data to a surface 
height field h( r, t) and to quantify properties of the deviations of the surface height field 
Ah(r,t) = h(r,t) — ho and surfactant density field Ar(r,t) = F(r,i) — r 0 . 

We perform all experiments just above the onset of the Faraday instability in order 
to excite a stable (non-chaotic) surface mode. To reach this regime, we quasi-statically 
increase the driving acceleration in steps of 1.5 x 10~ 3 g until the Faraday wave pattern 
is just identifiable in the moire images. Each quasi-static step during this preparation 
stage lasts for 90 s, during which the driving acceleration is held constant. In the data- 
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Figure 2. Stroboscopic imaging allows for the collection and phase-averaging of low-intensity 
fluorescence images over many cycles. During each second of the experiment, the container 
oscillates 20 times with an acceleration a(t) (black curve). During that second, the camera is 
triggered at four different phases of the oscillation (eight different phases of the Faraday wave 
oscillation). The eight distinct exposures (gray bars) are spaced 1 Faraday wave period apart 
with an exposure time of 1/8 of a Faraday wave period, illustrated by the width of the gray 
bars. The same stroboscopic pattern was used for both experiments, and both the FI and MI 
data. 


collection stage, the amplitude of the voltage signal for the electromagnetic shaker is 
held constant while we alternately collect MI and FI data. Although we collect data just 
above the Faraday wave onset, we nonetheless observe that the entire pattern rotates on 
the order of 10°/min. Similar rotations have been observed in other experiments, with 


the rotation axis slightly displaced from the center of the container (Gollub & Meyer 


1983 Douady et al. 1989). During image analysis, we remove the rotation by rotating 


images back to a common reference. 

To increase the signal-to-noise ratio, all fluorescence imaging is performed stroboscop- 
ically to phase-average over the low-intensity signal. The trigger pattern for the camera 
is chosen so that we cyclically sample eight different phases [i = 1... 8) of the Faraday 
waves (four phases of the meniscus waves), as shown in figure [2] In order to measure sub¬ 
period dynamics, we set the exposure time for each image to be 1/8 of the subharmonic 
Faraday oscillation period, and average 700 stroboscopically-acquired images at each of 
the eight phases. Obtaining this fluorescence composite data takes 14 minutes in total, 
during which time the Faraday wave rotates through roughly 140°. We take a series of 
32 moire images during each minute of the data collection to account for the rotation 
and to measure the surface height field. 


2.1. Apparatus and Materials 

In order to produce stable Faraday waves, the vertical driving requires small amplitude 
steps during the preparation stage, small drift in the acceleration amplitude during the 
data-collection stage, and pure sinusoidal driving at all times. We achieve these require¬ 
ments with a PM50A MB Dynamics shaker (peak force 220 N) mounted on a vibration 
isolation table. The driving acceleration is transferred from the shaker to the container 
by two parallel shafts. Each shaft is guided by a low-friction air bushing, and the parallel 
placement of the two bushings suppresses rotational motion. The mounting bar for the 
container contains 2 two-axis ADXL203 accelerometers, sampled at a rate of 20 kHz with 
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a noise threshold of approximately 10 _3 <?. We achieve an uncertainty in ao of 5 x 10~ 3 g 
and all higher harmonics are at least a factor of 10 -3 smaller than oq. 

The choice of materials in our experiments is guided by two important considerations: 
the sensitivity of surface experiments to preparation and handling, and the requirement 
of a low fluorescence background. The black-anodized container is made of aluminum 
machined into a cylindrical well of radius 14.6 cm. A 200 mm silicon wafer is embedded 
into the base of the well so that the wafer’s top surface is flush with the bottom of the 
well. The silicon wafer provides both a reproducible substrate and excellent reflectivity 
( Strickland et aL|2014 |. 

Prior to each experiment, we clean the container with detergent and perform a final 
rinse with 18.2 Mfi water before drying with dry nitrogen gas. The silicon wafer is 
cleaned for 5 minutes in an oxygen plasma environment, and all glassware is cleaned 
by soaking for several hours in a 2% Contrad 70 solution. These materials are rinsed 
with 18.2 Mil water and dried with dry nitrogen gas immediately before the experiment. 
To initialize each experiment, we fill the container to a depth of ho = 0.37 ± 0.02 cm 
(below the brim of the container) with 18.2 Mf2 water. Using a micropipette, we deposit 
a solution of chloroform and NBD-PC (l-palmitoyl-2-12-[(7-nitro-2-l,3-benzoxadiazol- 
4-yl)amino]dodecanoyl-sn-glycero-3-phosphocholine from Avanti Polar Lipids) onto the 
clean water surface. The concentration of NBD-PC in solution is 1 mg/mL, and due to 
the low interfacial tension between chloroform and water, the droplets spread over the 
fluid surface. 

In this paper, we report on two experimental runs: Experiment 1 with To = 0.3 /ig/cm 2 , 
and Experiment 2 with T 0 = 0.2/ig/cm 2 . Both runs are below the critical monolayer 
concentration for NBD-PC on water which is P c = 0.35/zg/cm 2 
The mean density To is calculated from the deposited volume of chloroform solution 
and the known dimensions of the container. To reduce disturbance from external air 
currents and dust, the entire apparatus is enclosed in a plastic tent. For Experiment 1, 
the temperature and humidity were 23.9 ± 0.4°C and 18.5 ± 0.5%, respectively; while for 
Experiment 2, they were 24.6 ± 0.4°C and 18.0 ± 0.5%. 


(Tsukanova et al. 2002). 


2.2. Fluorescence imaging 
The fluorescence imaging (FI) technique is based on quantifying spatiotemporal devia¬ 
tions of the surfactant density field (AT) from the mean density by observing the fluo¬ 
rescence intensity in digital images (|Fallest et al. 2010; Swanson et al. 2014 Strickland 


et al. 20l4| . The NBD-PC surfactant molecule contains a fluorophore with an excitation 


peak at 464 nm and an emission peak at 531 nm. Eight blue LEDs (1.5 W, 467 nm from 
Visual Communications Company, LLC) are mounted around the edge of the container 
to provide uniform excitation and illuminate the fluid surface at a low angle. The silicon 
wafer substrate reflects unabsorbed light away from the imaging system, thus minimiz¬ 
ing the background noise. An Andor Luca-R camera fitted with a Newport 530 ± 10 nm 
bandpass filter is positioned above the center of the experiment and takes images of sur¬ 
factant fluorescence within a field of view of width 15.4 cm. Because the camera images 
its own reflection near the center of the system, we report measurements for an annular 
region of the fluid surface. 

We calibrate the fluorescence intensity by performing experiments on a flat fluid surface 
covered with a known quantity of surfactant. In order to ensure a uniform density on 
the surface, we deposit the chloroform-dispersed surfactant onto the clean surface of 
water, wait 60 s for the chloroform to evaporate, and then drive the system at 20 Hz and 
0.2 g (above the Faraday onset) for 60 s. This promotes redistribution of the surfactant, 
so that after turning off the shaker and waiting for the fluid to settle back to its flat 
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Figure 3. Calibration of surfactant fluorescence intensity as a function of surfactant density. 
The stroboscopic exposure data were obtained with the same exposure time, number of images, 
and gain settings as used for the FI technique. The long exposure data was obtained with a 
single 2 s exposure and no gain, to confirm the trend. The markers show the mean intensity of 
the pixels within the field of view while the error bars illustrate the standard error of the pixels 
for the stroboscopic exposure data. The linear fit, shown as a dashed line, is used in the range 
r > 0.1 /rg/cm 2 of the experiments. 


state, the surfactant is uniformly distributed. In this state, we collect two types of data: 
700 images with the same exposure time as used for the FI composite images, and a 
final image with a single 2 s exposure time to provide a consistency-check on the trend. 
The mean and standard error of the fluorescence intensity as a function of T, reported 
in figure [3j are calculated from a composite of 700 images, as is done for each phase- 
averaged stroboscopically-acquired image. We observe that the calibration is linear for 
F > 0.1/xg/cm 2 . This corresponds to a regime in which the equation of state relating 
surfactant density and surface tension is also linear (jTsukanova et al. 20021. 

To obtain the FI composite images we average 700 individual stroboscopically-acquired 
images taken at the same temporal phase (see figure [2]), each one corrected for the 
accumulated pattern rotation. To measure AT, we first subtract a background image 
(obtained by averaging all eight phases) and then use the linear fit shown in figure [3] 
to convert the light intensity values to local surfactant density values at each pixel. A 
sample image is shown in figure [4] (a), with the brighter areas corresponding to higher 
surfactant density. 

2.3. Moire imaging 

In order to measure the surface height field, we use a moire imaging (MI) technique 
in which a reflected pattern of dark/light stripes reveals the spatial structure of the 
reflecting surface. For a known illumination pattern and known surface height field, a 
ray-tracing algorithm can produce a moire pseudo-image which should correspond to the 
observed image. To invert the moire image, we assume a physically-motivated functional 
form for the surface height field (see j |2.4[ ) and search for the optimal parameters in that 
ansatz which best reproduce all eight observed images. 

The moire imaging system consists of an Andor Luca-R camera which images the 
surface height field through a hole in the mask that creates the patterned light source (see 
figure [Tj) . The light source pattern is a transparent film printed with a set of concentric 
black rings of width 1.27 cm, chosen to be close to the wavelengths under consideration. 
This film, as well as a green gel filter (Lee Filters 736 Twickenham Green, with a peak 































Figure 4. (colour online) Sample FI composite image (a) showing AF from the i — 2 
stroboscopic phase of Experiment 1. sm oothed by a Gaussian filter of a = 0.04 cm. ( b ) The 
corresponding AIY, the best fit to \2.2\ . (c) The meniscus wave (harmonic) component and (e) 
the Faraday wave (subharmonic) component of the same Ar shown in (a). These components 
are isolated by adding and subtracting (and dividing by 2) the FI composite images from the 
i = 2 and i = 6 stroboscopic phases, (d) The meniscus wave (harmonic) component and (/) 
the Faraday wave (subharmonic) component of the same AFt shown in (6). ( g ) This image is 
generated by subtracting ( d ) from (c). The resulting pattern (smoothed with a Gau ssian filter 
of a = 0.09 cm) is a higher-order mode which is not accounted for in the ansatz (2.21. This 
mode has 24-fold symmetry and oscillates harmonically with peak-to-mean variation of 0.03 F c . 



Figure 5. Sample MI images, (a) Contrast-enhanced image and (6) binarization of the 
target pattern reflected off a flat fluid surface, (c) Contrast-enhanced moire image and ( d ) 
binarization of the target pattern reflected off the fluid surface at phase i — 2 for Experiment 
1. (e) Ray-traced pseudo-image and (/) binarization of the best-fit h{v) for image (c) using 
parameters given in table [T] 
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transmission wavelength of 525 ± 50 nm), are sandwiched between two panes of diffusing 
ground glass. Two 500 W incandescent bulbs positioned above the pattern illuminate the 
fluid surface through this target pattern, and the camera records the light reflected off the 
fluid surface. Even though the camera is focused on the fluid surface in order to optimize 
the low-intensity fluorescence measurements, the short focal length (35 mm) relative to 
the distance to the focal plane (67.8 cm) allows the use of a ray-tracing protocol to deduce 
the surface height profile. Sample images of the patterned light source, MI images, and 
best-fit ray trace are shown in figure [5j To reduce the number of rays that trace back to 
either the camera or outside the patterned light source, we use the same annular field of 
view considered for the FI images. 

Although the axes of the camera, pattern, and container are approximately aligned, we 
find it necessary to calculate their relative positions and orientations in order to achieve 
good ray-tracing results. To obtain the camera position and orientation, we image a set of 
metal posts of known height and spacing; to measure the center of the target pattern, we 
image the reflection off of the flat fluid surface. Using this data, we calculate the spatial 
resolution of the camera to be 0.0153 cm per pixel. To obtain the pixel coordinates of the 
system origin to sub-pixel resolution, we find the center of the meniscus wave pattern in 
composite data which has not been corrected for the Faraday wave rotation. The pattern 
in the rotation-uncorrected composite data is nearly identical to that in figure [4] (c). 

To find the best fit parameters of the ansatz for h(r,t), we maximize the cross¬ 
correlation between the MI images and ray-traced pseudo-images. The ray-tracing pro¬ 
tocol assumes that each ray corresponds to a pixel and all the rays start at a single 
point (the location of the camera) The rays are assumed to reflect off the fluid surface at 
the z = 0 plane, a justifiable assumption since the height of the wave peaks (O(10pm)) 
is much less than the distance between neighbouring peaks(0(lcm)). When waves are 
present, the normal to the surface changes orientation and the reflected ray responds 
accordingly. Each pixel of the ray-traced pseudo-images is assigned the intensity of the 
moire pattern at the pixel nearest the intersection of the ray and the plane of the pat¬ 
tern. For the few rays that trace back outside of the patterned light source, we assign 
the maximum intensity recordable by the camera. We account for motion blur due to 
the oscillation of the wave by averaging three ray-tracing results for each of the eight MI 
images. 


2.4. Surface height field and surfactant density field ansatz 

In order to compare Ah and AT, we use a physically-motivated ansatz which allows us 
to measure amplitudes, wavenumbers, and phases for the waves. In addition, we assume 
that Ah and AT are a linear superposition of the Faraday and the meniscus waves. Our 
results will indicate that this assumption accounts for the most significant signals in the 
data. To distinguish the ansatz from the measured MI and FI data, we use the subscript 
T for the ansatz quantities. 

Meniscus waves take the form of a superposition of inward and outward traveling 
Hankel functions (Saylor et al.\ 2000). The inward-propagating meniscus wave takes the 

form Re e lut (kr) where the complex-valued wavenumber k = + id measures 

both the spatial wavenumber kM and the damping factor a. Once at the center, the waves 
continue to propagate outward according to the form Re e lult H^\kr) . The resulting 
superposition is Re [e lujt Jo(fcr)]. 

Sinusoidally-driven Faraday waves on an uncontaminated surface of an infinitely deep 


inviscid fluid (Benjamin & Ursell 1954) or finite-depth low-viscosity fluid (Edwards & 


Fauve 


19941 take the form of a Bessel function J n (kpr) cos(nd) cos(^t) where u> is the 
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driving frequency and the spatial wavenumber k F is purely real. We allow for the Bessel¬ 
mode to be displaced from the center of the container, to account for the observed rota¬ 
tional instability, and use cylindrical coordinates measured with respect to the pattern- 
center rather than the axis of the experiment. 

Assuming a linear superposition of the meniscus and Faraday waves, our ansatz for 
Ah is given by 


A h T {r, 6 , t ) = h(r, 8 , t) - h 0 = A/i M t(u t) + A h FT {r', 8 ', t) 


= Ap[ Re 


e i("t+Mj oi (k H M + ia H )r) 


( 2 . 1 ) 




+ B h cos (—t + <f>HF)Jn(k HF r') cos (nd 1 


Sh) 


where the subscripts M, F indicate whether the parameter is associated with the meniscus 
wave or the Faraday wave, respectively. The primed coordinates are relative to the center 
of the Faraday wave instead of the axis of the system. 


Miles (1967) predicts that Ar F of a gravity-capillary wave is proportional to the 


Laplacian of the velocity potential, which is in turn proportional to A h F (Lamb 1945). 


Therefore, our AT ansatz takes the same form as (2.1): 


Ar T (r, 8 , t) = r(r, 8 , t)-T 0 = Ar„ T (r, t) + Ar FT (r', 8’, t) 


= ApRe 


e i(o,t+*™)j o(( fc rM+iar) r) 


( 2 . 2 ) 


+ B r cos (—t + 4‘rF)Jn(krFf') cos (nd 1 + <5 r ) 


Note that AT FT can be separated into independent spatial and temporal factors. In )j3j 
we will test whether the temporal dynamics of AT F follow the assumed cos(J (t) behavior. 


The nine fitting parameters in (2.1) (and the corresponding nine in (2.2)) each have 
a physical interpretation. To distinguish the parameters of A h F and ATp, we denote 
their parameters with the subscript H , T respectively. The amplitudes of the meniscus 
and Faraday waves are given by Ah and Bp, respectively. The spatial wavenumbers of 
the pattern are: k F M and kn f (real, radial for both waves), a (imaginary, radial for 
meniscus wave), n (real, azimuthal for Faraday wave). The value of a corresponds to the 
energy damping rate of the meniscus wave. Only the relative phase between AF and Ah 
is independent; we report temporal phase shifts for both the meniscus wave (4>hm — 
and the Faraday wave (</>hf — <t>rF)- Finally, there is also a spatial phase shift for the 
Faraday waves (Sh — hr). 

Due to the high-dimensionality of the parameter space, finding a best fit to this ansatz 
takes place in several stages, described below. The full set of fitting parameters for both 
A h F and AT F are provided in table [l] and discussed in detail in 

First, we decompose the FI images into a harmonic component (to fit the meniscus- 
wave terms) and a sub-harmonic component (Faraday terms) (see figure [4] (c,e)). This 
can be done by simply adding or subtracting images that are precisely one half Faraday 
wave period apart (e.g. images taken at phase i = 1 added to/subtracted from images 
taken at i = 5). We use a Levenberg-Marquardt algorithm to fit the resulting meniscus- 
series and Faraday-series to the corresponding terms in \2.2\ . A sample comparison of 
the FI composite data and the best-fit AT F are shown in figure [4] (6,d,/). 

Solving the inverse problem using the MI images is more challenging, and starts from 
an initial A h F chosen using the fluorescence-determined wavenumbers (with the missing 
parameters chosen by hand). Using MATLAB fminsearch , we simultaneously maximize 
the cross-correlation of all eight MI images against their corresponding ray-traced pseudo- 
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Experiment 1 

Experiment 2 


fluid thickness 

h 0 (cm) 

0.37 

0.37 


mean surfactant concentration 

To (f£g/cm 2 ) 

0.3 

0.2 


acceleration amplitude 

ao (g) 

0.145 

0.140 


symmetry number 

n 

12 

11 


meniscus: amplitude 

Ar (F c ) 

0.22 

0.04 

h 

meniscus: wavenumber 

krM (cm -1 ) 

6.70 

6.56 


meniscus: damping 

ar (cm -1 ) 

-0.15 

-0.30 


Faraday: amplitude 

B r (r c ) 

0.45 

0.20 


Faraday: wavenumber 

kvF (cm -1 ) 

3.40 

3.24 


meniscus: amplitude 

A h (cm) 

0.0022 

0.0033 

E-h 

meniscus: wavenumber 

k H M (cm -1 ) 

6.6 

6.5 


meniscus: damping 

a H (cm -1 ) 

-0.17 

-0.11 


Faraday: amplitude 

B h (cm) 

0.021 

0.021 


Faraday: wavenumber 

k HF (cm -1 ) 

3.4 

3.3 


meniscus: temporal phase shift 

4>HM — <t>TM 

1.19 

2.03 


Faraday: temporal phase shift 

4>HF — (j>TF 

2.38 

2.42 


Faraday: spatial phase shift 

8h — 5r 

0.33 

0.76 


Table 1. Summary of experimental parameters and best-fit parameters ((2.11 and (2.2)) for 
the two experiments. The uncertainty of all parameter values is in the last reported digit. The 
center of the Faraday wave (axis of patten rotation) is within ±3.5 mm of the container axis. 


images. To minimize the effects of pattern noise on the cross-correlation, we binarize and 
then Gaussian blur both the MI image and the ray-traced pseudo-image before perform¬ 
ing the cross-correlation. The searching algorithm proceeds by varying the parameters in 
) until a best fit is found. A sample image is shown in figure [5j illustrating that we 
are able to fit even small features of the surface waves. 


( 2 . 


3. Results 

We present direct measurements of the spatiotemporal dynamics of a surfactant mono- 
layer on gravity-capillary waves. In figures [4] and [5j we showed that our measurements 
of deviations of the surface height field Ah and surfactant density field AT are well- 
approximated by a linear superposition of a Faraday wave and a meniscus wave, each 
represented by a Bessel function. © and ( |2.2| ), together with the parameters provided 
by fits to the data (see table [l]), allow us to examine the spatiotemporal dynamics in 
detail. 

In figures [6] [8j we present three-dimensional visualizations of the complete Ah and AT 
fields (a snapshot at a single phase), as well as the same data split into its meniscus 
and Faraday wave components. In each case, the best-fit Ahp is shown as a surface 
mesh, with the colouring representing the measured AT: bright green regions are rich 
in surfactant, while dark green are depleted. Movie versions of these figures, showing 
all eight phases, are available in the online supplemental material. Below, we compare 
observations about both the traveling meniscus waves and the standing Faraday waves, 
along with the weaker higher-order mode observed in the FI composite data. 

For the meniscus waves, we observe that the surfactant accumulates on the leading 
edge of the traveling wave, as shown in figure[7](c), where the dashed line (ATmt) leads 
the solid line (AIimt)- The fields Ah and AT both exhibit a Bessel Jo mode, and the real 
parts of the associated wavenumbers k{jM and kpM agree, as expected. We also observe 
























Figure 6. (colour online) Visualization of A hr a nd Ar for Experiment 1, within the annular 
field of view. The best fit using the h-ansatz (2.11 (mesh surface) and the FI composite data 
(colouration) for the i = 2 stroboscopic phase. The peak-to-mean variations in A hr and AFt 
are ±60 /im and ±0.15 r c respectively. The AT fields are smoothed by a Gaussian filter of 
a = 0.04 cm to reduce noise. In figures [7] and [5] this same data is decomposed into into its 
meniscus and Faraday wave components. A video version of this plot is available in the online 
supplemental material. 


the expected sinusoidal evolution predicted by Saylor et al. (2000) and Picard & Davoust 


(2006). The maximum peak-to-mean variations in AhMT, measured within our annular 
field of view, are 5 and 8 /im for Experiments 1 and 2, respectively. The corresponding 
values for A T mt are 0.05 and 0.02 P c . These variations in A T mt correspond to variations 
in the surface tension of 0.5 and 1.5 mN/rn respectively (Tsukanova et al. 20021. The 
location of the surfactant relative to the fluid wave corresponds to a temporal phase shift 
between A/imt and ATmt, shown in figure [?] (c). 

For the Faraday waves, we observe that the surfactant accumulates in the troughs of the 
standing wave, as shown in figure[8](d), where the dashed line (AT ft) leads the solid line 
(A hpr)- The fields Ah and AT both exhibit a Bessel J n mode, and the wavenumbers knF 
and k rF agree, as expected (Miles 1967). We note that the wavenumber of the Faraday 
wave is roughly double that for the meniscus wave, and the values of the wavenumbers can 
be captured by the finite depth Kelvin dispersion relation ( uj 2 = (gk ± ak 3 /p) tanh(/iofc)) 
when the accepted values for the acceleration due to gravity g and density of water p 
are used and the same surface tension is left as a fitting parameter. The best fit surface 
tensions for Experiments 1 and 2 are 33.3 and 35.6 mN/rn respectively. The fields Ah and 
AT both exhibit a sinusoidal azimuthal behavior with the same symmetry n as well as a 
sinusoidal evolution at half the driving frequency. The maximum peak-to-mean variations 
in A hpTi measured within our annular field of view, are 62 and 64 pin for Experiments 
1 and 2, respectively. The corresponding values for ATft are 0.13 and 0.06 T c . These 
variations in AT ft correspond to variations in the surface tension of 1.5 and 3.5 mN/m 
respectively (Tsukanova et aZ.|;2002). Similar magnitudes in the variations of the surface 
tension have been observed for 1 Hz gravity waves contaminated with various molecular 
monolayers (Lange & Hiihnerfuss 1984). The location of the surfactant relative to the fluid 
wave corresponds to a temporal phase shift between Ah ft and ATp^r, shown in figure [8] 
( d ). We also observe a small but non-negligible spatial phase shift ( Sh — Sr) between 
AT ft and AhpT equal to 0.33 and 0.76 radians for Experiment 1 and 2 respectively. 
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AhMT " - AfMT 


FI Data 


Figure 7. (colour online) Menisc us w ave component of A hr and AF for Experiment 1. 
(a) The best fit AIimt ansatz from (mesh surface) and the harmonic component of AF 
(colouration) for the i = 2 stroboscopic phase. ( b ) Radial slice of both fields along d = 2.83 ra¬ 
dians, the orange strip shown in (a), at phase i = 2. The solid line is the best-fit height field 
A hMr(r), and the surfactant density field is shown both as the data AT m (or ange point s) and as 
the ansatz AFmt (dashed line). The lines are the Bessel modes given in (2.11 and ( |2.2[ ). (c) The 
temporal evolution of AhMT (solid line) and AYmt (dashed line) at the location r = 4.79 cm 
and 6 = 2 .83 r adians, marked by the magenta dot in (a). The lines are the sinusoids given in 
(2.11 and (2.2). In all cases, the AF fields are smoothed by a Gaussian filter of a = 0.04 cm to 
reduce noise. A video version of this plot is available in the online supplemental material. 


This spatial phase shift is in the same direction as the rotation of the Faraday wave 
pattern in both experiments. 

We have also considered a stricter test of the evolution for AT of the Faraday waves 
to determine whether they evolve non-sinusoidally, as indicated by the numerical solu¬ 


tions of Ubal et al. (2005). As given in (2.2), the ansatz for ATp can be factored into 


Sxir', 0') Tr{t) where the spatial factor is given by Sxir', 9') = Br </ ? ^(fc^F?' , ) cos (n9' + 
<5r) and the temporal factor by Tr{t) = cos(^f + (j>r f)- The subscript T denotes these 
quantities are fits to the ansatz, rather than the data itself. Since the spatial factor is 
well-characterised by the ansatz ( S = St), as shown in figures a (e,/) and [8] (&,c), we 
can isolate the temporal evolution of AT F for the data by considering the quantity 


m = 


AT F (r',9',t) 
S T {r', 6') 


(3.1) 


where the angle brackets denote an average over all space (approximately 5 x 10 5 pixels 
in total). In performing the average, we exclude pixels where \ST(r\0')\ < 0.006 r c to 
avoid the calculation being badly-conditioned. As shown in figure [8] (d), the average 
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Figure 8. (colour online) Far aday wave component of Ahr and Ar for Experiment 1. (a) 
The best fit /ift ansatz from (2.11 (mesh surface) and the sub-harmonic component of AT 
(colouration) for the i = 2 stroboscopic phase. (6) Radial slice of both fields along 9 C = 2.83 ra¬ 
dians, the orange strip shown in (a), at phase i = 2. The solid line is the best-fit height field 
AhFT{r,8 c ), and the surfactant density field is shown both as the data AIV (oran ge p oints ) 
and as the ansatz AVft (dashed line). The lines are the Bessel modes given in (2.1l and (2.2). 
(c) Azimuthal slice of both fields along r c = 4.26 cm, the cyan strip shown in (a), at phase 
i = 2. Here, the surfacta nt de nsity field data AF_f is represented with cyan points. The lines 
are the sinusoids given in ( |2. 1| ) and \2.2\ and plotted within a restricted range of the azimuthal 
dimension[— centered about t) c . ( d ) The temporal evolution of AIift (solid line) and 

AF_ft (dashed line) at the location r c cm and 6 C radians, marked by the magenta dot in (a). 


The lines are the s inusoids given in (2.1) and (2.2). The magenta circles in ( d ) plot behavior of 
T{t) given in ( |3. 1[ ), with errorbars much smaller than the symbol size. The AT fields in ( a,b,c ) 
are smoothed by a Gaussian filter of a = 0.04 cm to reduce noise. A video version of this plot 
is available in the online supplemental material. 


temporal evolution T{t) (and therefore the temporal evolution of ATp) is well-described 
by a sinusoidal function. 

Our ansatz assumes a linear superposition of the meniscus and Faraday waves. While 
this ansatz captures most of the total pattern, there is a residual mode in the FI data 
(see figure [4] ( g )) which generates peak-to-mean variations in AT of 0.03 and 0.01 r c 
for Experiments 1 and 2 respectively. This higher-order mode likely arises through a 
nonlinear interaction between the meniscus and Faraday waves. Although this mode 
has a Bessel-mode shape like the Faraday waves (F(r, 6) ~ J m (kr) cos {nd + </>)) with 
m = 2n (twice the symmetry number of the fundamental mode), the distance between 
neighbouring extrema is comparable to that of the meniscus waves. Additionally, the 
higher mode is harmonic in time, unlike the subharmonic Faraday wave. 

We observe that the degree of surfactant mobility depends on the type of wave. By 
comparing the T-amplitudes of the associated Faraday vs meniscus Bessel modes ( ^ ) > 
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we find a ratio of 2 for Experiment 1, and 5 for Experiment 2. This indicates that 
Faraday waves have a larger effect than meniscus waves. However, if the T-amplitudes are 
rescaled by the associated h- amplitudes, then meniscus have proportional larger effect: 


the ratio is 0.21 for Experiment 1 and 0.78 for Experiment 2. The cause of these 

differences is unclear and could be due to the difference between traveling vs. standing 
waves, to the difference in oscillation frequency, or to some subtlety of nonlinearity. In 
contrast, the observed differences between the two experiments is clear and indicates that 
material mobility depends on To- 

There are other subtle differences between the two experiments, suggesting an array 
of rich dynamics that depend on the mean surfactant concentration. First, the temporal 
phase shifts Af> = cf>H — 4>r are consistent for the Faraday waves across both experiments, 
but differ for the meniscus waves. Furthermore, only A (f>M for Experiment 1 is less 
than 7t/2, the expected maximum temporal phase shift for linear gravity-capillary waves 
(Lucassen-Reynders & Lucassen 1970). Second, the two damping factors (a# and or, 
both associated with the meniscus wave) agree for Experiment 1 but not Experiment 2. 
From the theoretical predictions of Picard & Davoust (2006), we would expect agreement. 


4. Conclusions 

We have successfully developed a novel technique for measuring the surface height h 
and surfactant density T fields for waves propagating on a surfactant-covered fluid. We 
observe that the surfactant accumulates on the leading edge of traveling meniscus waves 
and the troughs of standing Faraday waves. The deviations of both the surface height field 
Ah and surfactant density field Ar take characteristic Bessel forms for both the meniscus 
and Faraday waves. For the meniscus waves, both fields follow a Jo Bessel mode, but are 
temporally phase shifted from each other. For the Faraday waves, both fields follow the 
same J n Bessel mode, but are both temporally and spatially phase shifted from each 
other. 


For meniscus waves, Saylor et al. (2000) and Picard & Davoust (2006) have analytically 


derived Ah and AT. These two fields are expected to share the same functional form 
up to a complex valued coefficient (F(r, t) = F 0 Jo(fcr) e lut where k = kM + icx) and 
therefore are expected to have the same wavenumbers and damping factors. However, 
this work does not predict the relationship between the values of Jo for the two fields, 
which would be necessary to understand the temporal phase shifts. Lucassen- Reynders| 
& Lucassen ( |1970 ) studied 2-dimensional traveling gravity-capillary waves in a Cartesian 
geometry, and derived a prediction for the phase shift between the surface area expansion 
and Ah. For T inversely proportional to the area of the fluid surface, the temporal phase 
shift between Ah and AF is predicted to fall between 0 and 7 t/ 2 depending on the sur¬ 
face compression modulus. This corresponds to the surfactant accumulating somewhere 
between the crests of the waves and the leading edge. 

In Experiment 2, we observe </>#m — 4>tm > 7t/ 2 and an 7^ ar, both of which sug¬ 
gest that we may be exciting longitudinal waves in the system. Longitudinal waves, like 
gravity-capillary waves, are a solution to the linearized Navier-Stokes equations for in¬ 
compressible fluids. By including the effect of the surfactant through the normal and 


tangential stress boundary condition, Lucassen (1968a) derived a dispersion relation for 


which one branch corresponds to gravity-capillary waves, which have roughly equal parts 
transverse and tangential motion of the interface, and the other branch corresponds 
to longitudinal (i.e., Marangoni) waves, which have significantly more tangential mo¬ 
tion than transverse (Lucassen-Reynders & Lucassenjl970). The possibility of a gravity- 
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capillary wave resonantly exciting a longitudinal wave was first suggested by |Lucassen ; 
(1968a) because the maximum in energy dissipation coincides with the equality of the 
magnitude of the gravity-capillary and longitudinal wavenumbers. Recently, |Ermakov] 
(2003) proposed a mechanism for this resonant excitation which would explain the co¬ 
incidence of damping with the wavenumber equality. Future experiments which measure 
Ar would be able to test this mechanism. 

The discrepancies in the measured values of a highlight the importance of quantifying 
the surface dilational viscosity, dilational elasticity, and surface tension (collectively, the 
interfacial rheology). Historically, techniques for measuring the interfacial rheology have 
relied on extracting the wavenumber and damping factor from the surface height field h 
of a traveling gravity-capillary wave (Miyano et a(.|[l983 Jiang et aZ.||1993. Saylor et al. 


|2000 Behroozi et al. 2007). The experiments described here provide a way to make a 


more direct measurement, using both the h and T fields. 

For Faraday waves, there is theoretical literature relating the interaction between a 


surfactant monolayer and a Faraday wave in a fluid of arbitrary depth (Kumar & Matar 


2002 2004 Martin & Vega 2006 Ubal et al. 2005). Kumar & Matar (2002 20041 assumed 


that in an infinitely broad fluid, both Ah and AT oscillate sinusoidally and with zero 
temporal phase shift. They also predicted a spatial phase shift between the two fields. 
In the experimental data, we observe sinusoidal oscillations, but with a temporal phase 
shift of approximately 2.4 radians. However, we confirm the spatial phase shift (Sh — 
dr m table [l]). Martin & Vega (20061 attribute rotational drift of Faraday waves to 
the generation of a streaming flow, under the assumption that the two fields have no 
spatial or temporal phase shift. However, it is also plausible that our observed rotation 
of the Faraday waves is caused by the spatial phase shift. Ubal et al. (2005) studied the 
motion of Ah and AT through numerical simulations of the two dimensional Navier- 
Stokes equations with physically realistic interfaces and spatially symmetric boundary 
conditions in a finite depth fluid. In agreement with our experiments, they predicted 
that Ah and AT have similar spatial patterns, and that maxima in AT would precede 
maxima in Ah, an effect characterised by a temporal phase shift. The prediction of peak- 
to-mean variations in AT ss 0.03 to 0.13T C is consistent with our observations, although 
for different parameters than used in the experiment. However, their conclusion that 
the two fields have different non-sinusoidal time dependence does not agree with our 
data. Additionally, the use of spatially symmetric boundary conditions in the numerical 
simulations excludes the possibility of predicting the drift of the Faraday wave that we 
observe. 

Finally, although we have studied the wave-driven accumulation of a molecular layer 
of surface contamination, similar effects are known to be present for other types of sur¬ 
face contamination. For example, Sanl et al. (2014) used Faraday waves to redistribute 
a monolayer of millimeter-scale polystyrene spheres deposited on the surface. Depending 
on the concentration of the particulate layer, the macroscopic contaminant would accu¬ 
mulate in either the anti-nodes (peak/trough) or the nodes of the waves. This behavior 
is quite different from what we observe at the molecular scale, suggesting the presence 
of a crossover in particle size or surface activity. Because oceanic contamination ranges 
in size from the molecular scale to flotsam and jetsam, more work is needed to develop 
a complete picture of the dynamics of contamination on surface waves. 
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